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Abstract 

We study a restricted class of self-avoiding walks (SAW) which start at the origin (0, 0), end 
at {L,L), and are entirely contained in the square [0, L] x [0,L] on the square lattice Z^. The 
number of distinct walks is known to grow as 

+o{L )^ estimate A 1.744550 ± 0.000005 
as well as obtaining strict upper and lower bounds, 1.628 < A < 1.782. We give exact results for 
the number of SAW of length 2L + 2K for K — 0,1,2 and asymptotic results for K = o{L^/^). 

We also consider the model in which a weight or fugacity x is associated with each step of 
the walk. This gives rise to a canonical model of a phase transition. For x < the average 
length of a SAW grows as L, while for a; > 1/jU it grows as L^. Here fi is the growth constant 
of unconstrained SAW in Z^. For x = l//x we provide numerical evidence, but no proof, that 
the average walk length grows as L^/^. Another problem we study is that of SAW, as described 
above, that pass through the central vertex of the square. We estimate the proportion of such 
walks as a fraction of the total, and find it to be just below 80% of the total number of SAW. 
We also consider Hamiltonian walks under the same restriction. They are known to grow as 
+o{L ) ^j^p same L x L lattice. We give precise estimates for r as well as upper and lower 
bounds, and prove that r < A. 



1 Introduction 

We consider the problem of self-avoiding walks on the square lattice Z^. For walks on an infinite 
lattice, it is generally accepted that the number of such walks of length n, equivalent up to 
a translation, denoted 

grows as ~ const. fi^Ti^ , with metric properties, such as mean- 
square radius of gyration or mean-square end-to-end distance growing as {B?')n ~ const.n?'^ , where 
7 = 43/32 and = 3/4. The growth constant fi is lattice dependent, and for the square lattice is not 
known exactly, but is indistinguishable numerically from the unique positive root of the equation 
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13x^ — 7x'^ — 581 = 0. We denote the generating function by C{x) := X^^Cnx". It will be useful 
to define a second generating function for those SAW which start at the origin (0, 0) and end at a 
given point {u,v), as G'(o,o;n,?))(^)- terms of this generating function, the mass m{x) is defined 
[H] to be the rate of decay of G along a coordinate axis, 

I ^ 1- -logg(0,0;n,0)(^) 

mix) := lim ^ . (Ij 

n^oo n 

Here, we are interested in a restricted class of square lattice SAW which start at the origin (0, 0), end 
at (L, L) , and are entirely contained in the square [0, L] x [0, L] . A fugacity, or weight, x is associated 
with each step of the walk. Historically, this problem seems to have led two largely independent 
lives. One as a problem in combinatorics (in which case the fugacity has been implicitly set to 
X = 1), and one in the statistical mechanics literature where the behaviour as a function of fugacity 
X has been of considerable interest, as there is a fugacity dependent phase transition. 

The problem seems to have first been seriously studied as a mathematical problem by Abbott and 
Hanson ^ in 1978, many of whose results and methods are still powerful today. A key question 
considered both then and now, is the number of distinct SAW on the constrained lattice, and 
their growth as a function of the size of the lattice. Let Cn{L) denote the number of n-step SAW 
which start at the origin (0, 0), end at (L, L) and are entirely contained in the square [0, L] x [0, L\. 
Further, let Cl{x) := ^„c„(L)2;". Then Cl(1) is the number of distinct walks from the origin to 
the diagonally opposite corner of an L x L lattice. In pp, and independently in ^H], it was proved 
that Cl(1)^/^ A. The value of A is not known, though bounds and estimates have been given 
in ^^]. One of our purposes in this paper is to improve on both the bounds and the estimate. 

Like so many problems in lattice statistics, this one owes a debt to J. M. Hammersley. A closely 
related problem to the one considered here is discussed in ■, which is in turn devoted to problems 
posed by Hammersley. However the earliest mention of this problem appears to be by Knuth |12j . 
who calculated the number of SAW crossing a 10 x 10 square by Monte Carlo methods, and 
estimated the number to be (1.6 it 0.3) x 10^^. It is now known, see Table|^below, that the correct 
answer is 1.5687.. x 10^^. A related problem was studied by Edwards in He considered SAW 
starting at a point denoted the origin with end point a distance L from the origin, and no other 
points at distance L or greater. Let g{L) denote the number of such SAW. Then Edwards proved 
that YiraL^oo g{L)'^^/ ^ ^ exists and lies between 2.3 and 5.0. In our notation, Edwards has proved 
that 1.53 < A < 2.24. Edwards also proved that the same limit holds for SAW from the origin to 
the boundary of any convex, bounded subset of 1?. His numerical work led him to suggest that A 
is about 1.77. Our best estimate, given below, is 1.744550(5). 

The problem of Hamiltonian paths on an L x M rectangular grid, going from (0, 0) to (L, M) 
has also been considered previously. Earlier work is described in where Collins and Krompart 
also give generating functions for the number of such paths on grids with M = 1,2,3,4,5. In JU] 
Jacobsen and Kondev gave a field-theoretical estimate of the growth constant for Hamiltonian SAW 
on the square lattice, which must fill a square, as 1.472801 it 0.00001. 

In the statistical mechanics literature, the problem appears to have been introduced by Whittington 
and Guttmann JH] in 1990, who were particularly interested in the phase transition that takes place 
as one varies the fugacity associated with the walk length. All walks on lattices up to 6 x 6 were 
enumerated, and the estimate A = 1.756 ± 0.01 was given. At a critical value, Xc the average walk 
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length of a path on an L x L lattice changes from G(L) to 0(L^), where we define G(x) as follows: 
Let a{x) and b{x) be two functions of some variable x. We write that a{x) = G(6(x)) as x — > if 
there exist two positive constants ki and K2 such that, for x sufficiently close to xq, 

Ki b{x) < a{x) < K2 h{x). 

In ^H] the critical fugacity was proved to be at least l//i, its value was estimated numerically and 
was conjectured to be Xc = and in J3] the conjecture was proved by Madras. 

The problem was subsequently taken up by Burkhardt and Guim [2j, who extended the enu- 
merations given in ^Sj to 9 x 9 lattices, and used their data to give the improved estimate 
A = 1.743 ± 0.005. By considering SAW as the — > limit of the O(A^) model of magnetism, 
Burkhardt and Guim show that the conjecture Xc = 1//^ made in JHl on numerical grounds follows 
directly, though this is not a proof, unlike the subsequent result of Madras JH] • 

They also gave a scaling Ansatz for the behaviour of Cl{x) for L large in the vicinity of x = Xc- 
They proposed 

CL{x)r.L-^^f[L'/''{x,-x)] (2) 



where v = 3/4, as described above, and r]c = 5/2 is the corner exponent of the magnetisation 
given by Cardy's [HI result r]c{9) = f for a wedge-angle 9, which is 7r/2 in this case. 7?|| = 5/4 is 
the surface exponent that characterises the decay of spin-spin correlations parallel to the boundary 
in the semi-infinite geometry, corresponding to wedge-angle vr. Consequences of this scaling Ansatz 
include the following predictions: 

Cl{xc) ~ const. L'"^" 
d 

{n{xc, L)) = x—Cl{xc) ~ const.L^/" (3) 

((n(x„L) - (n(x„L)))2) = (x^)^ hi Cl(x,) ~ const. . 

ox 

They tested these results from their numerical data, and found them well supported. We provide 
even firmer support for these results on the basis of radically extended numerical data. Equation 
(jnj has also previously been given by Duplantier and Saleur 0. 

Burkhardt and Guim also considered a generalisation of the problem considered here by including 
a second fugacity, associated with steps in the boundary. This allows the problem of adsorbing 
boundaries to be studied. We will not discuss this aspect of the problem further, except to note 
that in 12] a full scaling theory is developed, and the predictions of the theory are tested against 
numerical data. 



In ^ the slightly more general problem of SAW constrained to an L x M lattice was considered, 
where the analogous question was asked: how many non-self-intersecting paths are there from (0, 0) 
to (L, M)? If one denotes the number of such paths by Cl,m, it is clear that, for M finite, the paths 
can be generated by a finite dimensional transfer matrix, and hence that the generating function 
is rational ^7j- Indeed, in P it was proved that 

0,M = ^C.,.^ ^ (4) 

L>0 
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Figure 1: An example of a SAW configuration crossing a square (left panel), traversing a square 
from left to right (middle panel) and a cow-patch (right panel). 

(where here we have corrected a typographical error). It follows that Cl,2 ~ const. A2^, where 
A2 = y^= 1.81735.... 

In this paper we also consider two further problems which can be seen as generalisations of the 
stated problem. Firstly, we consider the problem where SAWs are allowed to start anywhere on the 
left edge of the square and terminate anywhere on the right edge; so these are walks traversing the 
square from left to right. We call such walks transverse walks. Secondly, we consider the problem 
in which there may be several independent SAW, each SAW starting and ending on the perimeter 
of the square. The SAW are not allowed to take steps along the edges of the perimeter. Such walks 
partition the square into distinct regions and by colouring the regions alternately black and white 
we get a cow-patch pattern. Each problem is illustrated in Figure ^ 

Following the work in ^HI? Madras ^H] proved a number of theorems. In fact, most of Madras's 
results were proved for the more general d-dimensional hyper-cubic lattice, but here we will quote 
them in the more restricted two-dimensional setting. 

Theorem 1 The following limits, 

fii{x) := lim Cl{x)^^^ and H2{x) := lim Cl{x)^^^\ 
are well-defined in M U {+00} . 
More precisely, 

(i) Hi{x) is finite for < x < l//.f, and is infinite for x > 1/fi. Moreover, < /ii(x) < 1 for 
< X < 1//U and ni^l/fi) = 1. 

(ii) fJ^2ix) is finite for all x > 0. Moreover, H2ix) = 1 for < x < l//.i and f^2ix) > 1 for x > 

In JH] the existence of the limit ii2{x) was proved, and in addition upper and lower bounds on 
H2{x) were established. 
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The average length of a (weighted) walk is defined to be 

(n(x,L)) := j;nc„(L)x"7^c„(L)x". (5) 

n n 

Theorem 2 For < x < we have that {n{x,L)) = G(L) as L ^ oo, while for x > l/fi, we 
have {n{x,L)) = 6(^2). 

In ^Hl it was proved that (n(l))L = Q{L'^). The situation at x = is unknown. We provide 
compelling numerical evidence that in fact {n{l/iJ,))L = @{L^^'^) , where u = 3/4, in accordance 
with an intuitive suggestion in both [2] and |13j . 

Theorem 3 For x > 0, define fi{x) = logiii{x) and f2{x) = logfi2{x). 

(i) The function f\ is a strictly increasing, negative-valued convex function of log x for < x < 
and fi{x) = 0(— m(x)) as x — > 1//^", where m{x) is the mass, defined by 

(ii) The function f2 is a strictly increasing, convex function of log x for x > l/^u, and satisfies 
< f2{x) < log fi + log X. 

Some, but not all of the above results were previously proved in JH]) but these three theorems 
elegantly capture all that is rigorously known. 

The rest of the paper is organised as follows: In the next section we describe our enumeration 
methods, and explain how they are used to obtain radically extended series expansions for the 
number of walks crossing a square, the number of cow-patch configurations and the number of 
transverse SAW. Section 3 details the results we have obtained. In Section 4 we derive methods 
for obtaining rigorous upper and lower bounds on A. In that section we show that upper bounds 
based on counting cow-patch configurations are fully equivalent to the method of Abbott and 
Hanson, based on 0-1 admissible matrices. An improved method of lower bounds based on counting 
transverse walks is also derived. In section 5 we then apply these methods to our radically extended 
enumerations to provide significantly improved bounds on A. In section 6 we give exact results for 
short SAW crossing a square. The shortest SAW that can cross a square from (0, 0) to (L, L) 
is of length 2L. We give the exact number of such SAW of length 2L + 2K, for K = 0,1, 2, and 
asymptotic results for K = o{L^^^). Section 7 is devoted to a numerical analysis which gives precise 
(though non-rigorous) estimates of A, for all three types of configurations, a discussion of the mean 
number of steps as a function of fugacity, fiuctuations in this quantity, and a scaling theory for such 
fluctuations. We also speculate on the nature of the sub-dominant behaviour of the asymptotic 
form for the number of SAW. Section 8 is also a numerical study, but of the number of SAW that 
pass through the central vertex of an L x L square. Finally in section 9 we study Hamiltonian 
paths, obtaining both rigorous upper and lower bounds on the growth constant, and a numerical 
estimate. 
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2 Exact enumeration 



In the following we give a fairly detailed description of the algorithm we use to enumerate the 
number of walks crossing a square and briefly outline how this basic algorithm is modified in order 
to include a step fugacity, study SAWs traversing a square and the cow-patch configurations. 

2.1 The basic algorithm 

We use a transfer matrix algorithm to count the number of walks crossing L x M rectangles. The 
algorithm is based on the method of Conway et al. |Sj for enumerating ordinary self-avoiding walks. 
The transfer matrix technique involves drawing a boundary line through the rectangle intersecting 
M + 1 or M + 2 edges. 

For each configuration of occupied or empty edges we maintain a count of partially completed 
walks intersecting the boundary in that pattern. Walks in rectangles are counted by moving the 
boundary adding one vertex at a time (see Figure [21) ■ Rectangles are built up column by column 
with each column constructed one vertex at a time. Configurations are represented by lists of states 
{cjj}, where the value of the state at position i must indicate if the edge is occupied or empty. 
An empty edge is indicated by ai = 0. An occupied edge is either free (not connected to other 
edges) or connected to exactly one other edge via a path to the left of the boundary. We indicate 
this by (Tj = 1 for a free end, cjj = 2 for the lower end of a loop and cr^ = 3 for the upper end 
of a loop connecting two edges. Since we are studying self-avoiding walks on a two-dimensional 
lattice the compact encoding given above uniquely specifies which ends are paired. Read from the 
bottom the configuration along the intersection in Figure |2l is {2203301203} (prior to the move) 
and {2300001203} (after the move). 

There are major restrictions on the possible configurations and their updating rules. Firstly, since 
the walk has to cross the rectangle there is exactly one free end in any configuration. Secondly, all 
remaining occupied edges are connected by a path to the left of the intersection and we cannot close 
a loop. It is therefore clear that the total number of 2's equals the total number of 3's. Furthermore, 
as we look through the configuration from the bottom the number of 2's is never smaller than the 
number of 3's (they are perfectly balanced parentheses). We also have to ensure that the graphs 
we construct have only one connected component. In the following we shall briefly show how this 
is achieved. 

We call the configuration before and after the move the 'source' and 'target', respectively. Initially 
we have just one configuration with a single '1' at position (all other entries '0') thus ensuring 
that we start in the bottom-left corner. As the boundary line is moved one step, we run through 
all the existing sources. Each source gives rise to one or two targets and the count of the source is 
added to the count of the target (the initial count of a target being zero). After a source has been 
processed it can be discarded since it will make no further contribution. Tabled lists the possible 
local 'input' states and the 'output' states which arise as the kink in the boundary is propagated 
one step and the various symbols will be explained below. 

Firstly, the values of the 'Bottom' and 'Top' table entries refer to the edge-states of the kink prior 
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Figure 2: The left panel shows a snapshot of the intersection (dashed line) during the transfer 
matrix calculation. Walks are enumerated by successive moves of the kink in the boundary, as 
exemplified by the position given by the dotted line, so that the L x M rectangle is built up one 
vertex at a time. To the left of the boundary we have drawn an example of a partially completed 
walk. Numbers along the boundary indicate the encoding of this particular configuration. The 
right panel shows some of the local configurations which occur as the kink in the intersection is 
moved one step. 



Table 1: The various 'input' states and the 'output' states which arise as the boundary line is 
moved in order to include one more vertex. Each panel contains up to three possible 'output' states 
or other allowed actions. 
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to the move. The Top (Bottom) entry is the state of the edge intersected by the horizontal (lower 
vertical) part of the boundary. 

Some of the updating rules are illustrated further in Figure [2 The topmost panels represent the 
input state '00' having the allowed output states '00' and '23' corresponding to leaving the edges 
empty or inserting a new loop, respectively. The middle panels represents the input state '20' with 
output states '20' and '02' from the two ways of continuing the loop end (note that the loop has to 
be continued since we would otherwise generate an additional free end not located at the allowed 
positions in the corners). The bottommost panels represents the input state '22' as part of the 
configuration {02233}. In this case we connect two loop ends and we thus join two separate loops 
into a single larger loop. The matching upper end of the innermost loop becomes the new lower 
end of the joined loop. The relabeling of the matching loop-end when connecting two '2's (or two 
'3's) is denoted by over-lining in Tabled When we join loop ends to a free end (inputs '12', '21', 
'13', and '31') we have to relabel the matching loop end as a free end. This type of relabeling is 
indicated by the symbol 00. The input state '11' never occurs since there is only one free end. 
The input state '23' is not allowed since connecting the two ends results in a closed loop and we 
thus discard any configuration in which a closed loop is formed. It is quite easy to avoid forming 
closed loops. We only have to be careful when the input is '03' or '30'. If the upper end of the 
loop is continued along the vertical output edge we would form a closed loop if the horizontal edge 
immediately below was a lower loop end, and we just check the state of this edge and only proceed 
if it is not in state '2' (naturally the upper loop end can always by continued along the horizontal 
output edge). 

Finally, we have marked two outputs, from the inputs '01' and '10' with 'Res', indicating situations 
where we terminate free ends. This results in completed partial walks and is only allowed if there 
are no other occupied edges in the source (otherwise we would produce graphs with separate pieces) 
and if we are at the top- most vertex (otherwise we would not cross the rectangle) . The count for 
this configuration is the number of walks crossing a rectangle of height M and length L equal to 
the number of completed columns. 

The time required to obtain the number of walks on L x M rectangles grows exponentially with 
M and linearly with L. Time and memory requirements are basically proportional to the maximal 
number of distinct configurations along the boundary line. When there is no kink in the intersection 
(a column has just been completed) we can calculate this number, A'conf(-^)> exactly. Obviously 
the free end cuts the boundary line configuration into two separate pieces. Each of these pieces 
consists of 'O's and an equal number of '2's and '3's with the latter forming a perfectly balanced 
parenthesis system. 

Each piece thus corresponds to a Motzkin path ^3 Ch. 6] (just map to a horizontal step, 2 to a 
north-east step, and 3 to a south-east step). The number of Motzkin paths M„ with n steps is easily 
derived from the generating function M{x) = ^„ MnX^, which satisfies M. = 1 + xM. + x'^M.'^ , so 
that 

M{x) = [l-x- V(l + a;)(l - 2,x)]/2x^. (6) 

The number of configurations A'^conf (-'^) for a rectangle of height M is simply obtained by inserting 
a free end between two Motzkin paths, so that the generating function Y^j^^ Nconii^)^^ is simply 
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xM{x)'^. The La grange inversion formula gives 



When the boundary Une has a kink the number of configurations exceeds Ncon{{M) but clearly is 
less than A'conf (-^+ !)• From © we see that asymptotically Nconi{M) grows like S*'^ (up to a power 
of M). So the same is true for the maximal number of boundary line configurations and hence 
for the computational complexity of the algorithm. Note that the total number a walks grows like 
, so our algorithm leads to a better than exponential improvement over direct enumeration. 

The integers occurring in the expansion become very large so the calculation was performed using 
modular arithmetic This involves performing the calculation modulo various prime numbers 
Pi and then reconstructing the full integer coefficients at the end. We used primes of the form 
Pi = 2^^ — Vi where rj are distinct integers, less than 1000, such that pi is a (different) prime for 
each value of i. The Chinese remainder theorem ensures that any integer has a unique representation 
in terms of residues. If the largest integer occurring in the final expansion is m, then we have to 
use a number of primes k such that piP2 • • • Pfe > m. 



2.2 Extensions of the algorithm 

The algorithm is easily generalised to include a step fugacity x. The count associated with the 
boundary line configuration has to be replaced by a generating function for partial walks. Since we 
only use this generalisation to study walks crossing an L x L square the generating function is just 
a polynomial of degree (at most) L{L + 2) in x. The coefficient of is just the number of partial 
walks of length n intersecting the boundary line in the pattern specified by the configuration. The 
generating function of the source is multiplied by x™" and added to the target, where m is the 
number of additional steps inserted. Not all L{L + 2) terms in the polynomials need be retained. 
Firstly, any walk crossing the square has even length. Thus in the generating functions for partial 
walks either all the even or all the odd terms are zero, and we need only retain the non-zero terms. 
Secondly, in order to construct a given boundary line configuration, a certain minimal number of 
steps njnin are required. Terms in the generating function of degree lower than rimin are therefore 
zero and again we need not store these. 

The generalisation to traversing walks is also quite simple. Firstly, we have M + 1 initial configura- 
tions which are empty except for a free end at position < j < M. This corresponds to the M + 1 
possible starting positions for the walk on the left boundary. Secondly, we have to change how we 
produce the final counts. The easiest way to ensure that a walk spans the rectangle and that only 
single component graphs are counted is as follows: When column L + 1 has been completed we 
look at the M + 1 configurations with a single free end and add the counts from all of them. This 
is the number of walks traversing an L x M rectangle. 

The generalisation to cow-patch patterns is more complicated. Graphs can now have many separate 
components each of which is a SAW, and there can thus be many free ends in a boundary line 
configuration. Note that each SAW starts and terminates with a step perpendicular to the border 
of the rectangle and there are never any steps along the edges of the borders of the rectangle. There 
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are 2 initial configurations since any of the edges in the first column from position 1 to M — 1 
can be occupied by a free end or be empty (recall that in cow-patch configurations the top and 
bottom-most horizontal edges cannot be occupied). There is an extra updating rule in the bulk in 
that we can have the local input '11' (joining of two free ends) with the only possible output being 
'00'. Also the updating rules at the upper and lower borders of the rectangle are different in this 
case. At the upper border we only have the input '00' with the outputs '00' and '10' corresponding 
to the insertion of a free end on a vertical edge at the upper border. There is no '23' or '01' outputs 
since these would produce an occupied edge along the upper border. At the lower border we have 
inputs '00', '01', and '02' and in each case the only possible output is '00' (with the appropriate 
relabeling in the '02' case). Finally, the count of the number of cow-patch patterns is obtained by 
summing over all boundary line configurations after the completion of a column. 



3 Results 



As discussed above, in order to obtain the exact value of the number of SAW crossing a square, 
some of which are integers with nearly 100 digits, we performed the enumerations several times, 
each time modulo a different prime. The enumerations were then reconstructed using the Chinese 
Remainder Theorem. Each run for a 19 x 19 lattice took about 72 hours using 8 processors of a 
multiprocessor 1 GHz Compaq Alpha computer. Ten such runs were needed to uniquely specify 
the resultant numbers. 

Proceeding as above, we have calculated c„(L) for all n for L < 17. In other words, we have obtained 
the polynomials Cl{x) for L < 17. In addition, we have computed Ci8(l) and Ci9(l), the total 
number of SAW crossing an 18 x 18 and 19 x 19 square respectively. We have also computed the 
corresponding quantities for cow-patch and transverse SAWs, denoted -Pl(I) and Tl(1) respectively, 
for L < 19. These are given in Tabled 

1 

In ^ the question was asked whether C^^/ decreasing in both L and M? We can answer this 
in the negative, based on our enumerations. 



4 Proofs of bounds 



Let C{L) be the set of self-avoiding walks crossing the Lx L square from its south-west corner (0, 0) 
to its north-east corner (L,L). Let C{L) denote the cardinality of C(L). Let T(L) be the set of 
self-avoiding walks that traverse, the L x L square: by this, we mean that the walk starts from the 
west edge of the square and ends on the east edge (Figure Let T(L) be the cardinality of T(L). 
Finally, let V{L) be the set of cow-patches, of size L: a cow-patch is a configuration of mutually 
avoiding self-avoiding walks on the L x L square, such that each walk has both endpoints on the 
border of the square, but never contains an edge of the border (Figure^. Let P{L) be the number 
of cow-patches of size L. 
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Table 2: The total number of walks crossing a square, Cl{1), cow-patch walks, -Pl(I) and traversing 
walks, Jl(1). 

L Cl{1) 

1 2 

2 12 

3 184 

4 8512 

5 1262816 

6 575780564 

7 789360053252 

8 3266598486981642 

9 41044208702632496804 

10 1568758030464750013214100 

11 182413291514248049241470885236 

12 64528039343270018963357185158482118 

13 69450664761521361664274701548907358996488 

14 227449714676812739631826459327989863387613323440 

15 2266745568862672746374567396713098934866324885408319028 

16 68745445609149931587631563132489232824587945968099457285419306 

17 6344814611237963971310297540795524400449443986866480693646369387855336 

18 1782112840842065129893384946652325275167838065704767655931452474605826692782532 

19 1523344971704879993080742810319229690899454255323294555776029866737355060592877569255844 

1 1 

2 7 

3 160 

4 11408 

5 2522191 

6 1718769373 

7 3598611604598 

8 23098353998190640 

9 453839082673896579243 

10 27266319759961440667165921 

11 5005013940387988257218110301496 

12 2805250606288167736619664411164848668 

13 4798636658841347169993094278185741344065154 

14 250425637137809429696661 10695844976426050692260400 

15 398585071868378544875200967972920693215965420927547891443 

16 19343509060397504009184634223201418820841655935064055180184148711 

17 2861743739297615012905209591294651941414000218185488280077237678797763881 

18 1290420684731131093964422300362403673911432011198730662653676329480448243238167005 

19 1773260101104126884305729846781529391070539884533101171392023893295633931250883380602647575 

1 8 

2 95 

3 2320 

4 154259 

5 30549774 

6 17777600753 

7 30283708455564 

8 152480475641255213 

9 2287842813828061810244 

10 102744826737618542833764649 

1 1 13848270995235582268846758977770 

12 5613766870113075134552249300590982081 

13 6856324633418315229580098999727214234534626 

14 25264653780547704599613926971040640439380254497299 

15 281194924965510769640501069703642937039678809002355743600 

16 9461739046646537749639494171503923182753987897972167546351180871 

17 963236702020101408274810653629921860636656580683490560257709270360444788 

18 296872411379358777499142156584947972393781613934413706389772635139720532797697401 

19 277150300263332125727926989254635730407844207233646123561354535935393720183262709640734296 
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L L-1 

Figure 3: From transverse walks to walks crossing a square (left). From walks crossing a square 
to cow-patches (right). 



We first prove in this section that 

limC(L)i/^' = limT(L)i/^' = limP(L)^/^' = A. (7) 
Then, we prove the following bounds on A: for L > 1, 

C(L)V{^+i)' < A, ^(^^i/((L+i)(L+2)) < \<{2P{L)f"^\ 

Let us first focus on Q. As recalled in the previous sections, the convergence of C(L)^/^^ to A has 
been proved in earlier papers ^|^. For walks of T{L), a similar result follows from the fact that 

C(L) < T{L) < C{L + 2). 

The first inequality above is obvious. The second one is explained on the left of Figure 01 

For cow-patches, the existence and value of limP(L)^/^^ follows from 

C{L-1) < P(L) < C{L + 3). 

The first inequality is explained on the right of Figure 01 The second one is a bit more tricky. We 
borrow the following argument from [7]. It is illustrated in Figurejll Start from a cow-patch of size 
L. Colour all cells of the square in black and white, in such a way that the south-west corner of the 
square is black and each step included in one of the walks of the cow-patch is adjacent to a black 
cell and a white one. Surround the square by a layer of black cells, so as to obtain a square of size 
L + 2, containing a certain number of white regions. For each white region, dig a tunnel (exactly 
one tunnel) in the outer layer to connect it to the outer world. In the figure thus obtained, the 
border of the black region forms a self-avoiding polygon, that includes each walk of the cow-patch. 
It remains to extend this polygon in a canonical way to obtain a walk of C(L -|- 3), illustrated in 
the last panel of Figure jll 

Let us now discuss lower and upper bounds on A. The left-hand side of Figure shows that for 
all £ and all odd k, it is possible to combine k'^ elements of C{£) to form an element of C{L) with 
L = k{l + 1). In Figure 01 k = 3. This shows that 

C(£)'=' < C{L). 
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L L 
Figure 4: From cow-patches to walks crossing a square. 



Hence 

C7(^)V{m)^ <C(L)i/^'. 
Taking the limit as k ^ co implies that for all £, 

< A. 

Similarly, let us try to pack transverse walks densely. The right-hand side of Figure shows that 
for all £ and /c, it is possible to combine A;^(£ + + 2) elements of T{£) to form an element of 
C{L) with L = k{e + l){£ + 2). This shows that 

T(£f('+^)('+^) < C{L). 

Hence 

2^(^)i/{(^+i){<?+2)) ^ c{Lf/^^. 
Taking the limit as A; — > oo implies that for all I, 

x>T{£y/mm+2))_ (8) 



Let us finally give upper bounds for A. Define a coloured cow-patch as a cow-patch in which the 
various regions are coloured in black and white, in such a way that two adjacent regions have 
different colours. Clearly, each cow-patch gives rise to 2 coloured cow-patches. Observe that there 
is a bijection between coloured cow-patches of size L and the admissible, matrices of the same size, 
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k{e+l) blocks 




as defined in SectionlHl Since an element of C{L), with L = ki, can be seen as tlie juxtaposition of 
k'^ admissible matrices (or colom'ed cow-patches) of size £, 

C{L) < {2P{l)f. 

That is, 

and by letting k oo, we obtain Abbott and Hanson's bound: for all £, 

One possible attempt to improve this bound is to consider generalised cow-patches, in which the 
walks are allowed to include edges lying on the west and south borders of the square (Figure inj. 
Let GP{L) denote the number of generalised cow-patches of size L. Since an element of C(L), with 
L = ki, can be seen as the juxtaposition of k'^ generalised patches, the above argument gives 

A < GP{lY'^\ 

We have not exploited this improvement, as it only changes the fourth significant digit of our bound. 



5 Bounds on the growth constant A 



For the more general problem of SAW going from (0, 0) to (L, M) on an L x M lattice, it was 
proved in P that 
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Figure 6: A generalized cow-patch. 



Theorem 4 For each fixed M , hnii^oo C^m ~ -^Af exists. 

1 

Further, Abbott and Hanson state that a similar proof can be used to establish that linii^oo C'/'^ •= 
A exists. This was proved rather differently in |18j . 

5.1 Upper bounds on A 

In ^ an upper bound on the growth constant A was obtained by recasting the problem in a matrix 
setting. We give below an alternative method for establishing upper bounds, based on defining a 
superset of paths. We then show that these two methods are in fact identical. 

Following consider any non-intersecting path crossing the L x L square. Label each unit 
square in the L x L lattice by 1 if it lies to the right of the path, and by if it lies to the left. 
This provides a one-to-one correspondence between paths and a subset of L x L matrices with 
elements or 1. Matrices corresponding to allowed paths are called admissible, otherwise they 
are inadmissible. Since the total number of L x L — 1 matrices is 2^^ , we immediately have the 
weak bound Cl,l < 2^ . Of the 16 possible 2x2 matrices, only 14 can correspond to portions of 
non-intersecting lattice paths. Note that there are only 12 actual paths from (0,0) to (2,2), but 
a further two matrices may correspond to paths that are embedded in a larger lattice. Thus we 
find the bound Cl,l < 14(^/^)2^ so A < 1.9343... Similarly, for 3 x 3 lattices we find 320 admissible 
matrices (out of a possible 512), so A < 320^/^ = 1.8982.. For 4x4 lattices, jTj claims that there 
are 22662 admissible matrices, but we believe the correct number to be 22816, giving the bound 
A < 1.8723... We have made dramatic extensions of this work, using a combination of finite-lattice 
methods and transfer matrices, as described below, and have determined the number of admissible 
matrices up to 19 x 19. There are 3.5465202 . . . x lO^'' such matrices, giving the bound 

A < 1.7817. 

This bound is fully equivalent to the bound A < {2Pi,)^/^^ , where Pl denotes the number of 
cow-patch configurations on the L x L lattice. This bound is proved below, in Section [IJ and the 
equivalence follows upon colouring cow-patches by two colours, such that adjacent regions have 
different colours. Labeling the two colours and 1 produces a — 1 matrix representation. 
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5.2 Lower bounds on A 



In P the useful bound 

M 

A > A-+^ (9) 

is proved. 

The above evaluation of A2, see @, immediately yields A > 1.4892 .... 

Based on exact enumeration, we have found the exact generating functions Gm{z) = '^Zl^lM'^^ 
for M < 6. For M = 3 we find: 

^ [1,-4,-4,36,-39,-26,50,6,-15,1] 

[1, -12, 54, -124, 133, 16, -175, 94, 69, -40, -12, 4, 1] ' 

where we denote by [oq, oi, . . . , a„] the polynomial + aiz + • • • + anz"^ . As explained above, all 
the generating functions Gm{z) are rational. For M = 4,5,6, their numerator and denominators 
are found to have degree (26,27), (71,75) and (186, 186) respectively, in an obvious notation. 

From these, we find the following values: A3 = 1.76331 . . ., A4 = 1.75146 . . ., A5 = 1.74875 . . . and 
Ag = 1.74728 .... Then from eqn. @ and Ag we obtain the bound A > 1.61339 .... 

However, an alternative lower bound can be obtained from transverse SAWs, defined in Section Q 
If Tl denotes the number of transverse SAW on the Lx L lattice, then we prove in the next section 
that 



A>r(L)i/«^+^)(^+2)). (10) 
From our enumerations of T(L), given above for L < 19, we obtain the improved bound A > 1.6284. 

Combining our results for lower and upper bounds finally gives 

1.6284 < A < 1.7817. 



6 Short walks crossing a square 



As defined in the introduction, let c„(L) be the number of n-step self-avoiding walks crossing an 
L X L square. Clearly, this number is zero when n is odd and also when n < 2L. It is almost as 
clear that 

.ML) ^ {f 

Indeed, there are 2L steps in the path, of which L must go north and L must go east. Note that 
the number C2l{L) has asymptotic expansion 

4^/1 1 

1 - + TTTTTTTT + .0 + 



Lvr V 4L 128L2 1024L3 
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Let us now prove that 

A walk counted by C2l+2(-^) has either L + 2 vertical steps (and L horizontal ones), or L vertical 
steps (and L + 2 horizontal ones). By symmetry, we can focus on the first case. Let w be such 
a walk. We say that lu has a vertical defect.. Among the L + 2 vertical steps of exactly one 
goes south, while the L + 1 others go north. The unique south step S is necessarily preceded and 
followed by an east step, which we denote respectively Ei and E2. Let us mark Ei and delete S 
and E2 (Figure [7j). The marked path w' thus obtained allows one to recover the original path w. It 
contains L + 1 north steps and L—\ east steps, one of which is marked. Moreover, the marked step 
cannot be at ordinate 0, nor at ordinate L + 1. Conversely, any walk w' satisfying these properties 
is obtained (exactly once) from a walk counted by C2l+2{L) and having a vertical defect. 




Figure 7: Enumeration of self-avoiding walks with one vertical defect. 

Marking one of the 

east steps gives a factor (L — 1). Now we must subtract the number of walks in which the marked 
step is either at level or at level + 1. Transforming the marked step into a vertical step shows 
that each of these two families of marked walks is in bijection with walks formed with L + 2 up 
steps and L — 2 down steps. Putting these observations together gives 

w.,^<(.-.,(/_\)-<-,))^..(-, 

Note that the number C2l+2{L) has the asymptotic expansion 

L4^ / 33 1345 23835 

2- — + —^-——r + 



\ 4L 64L2 512L3 



The same ideas may be used to find the value of C2L+4iL). We will prove that 

lc2L+4{L) = (48 + 90L + 8L2 - 28^^ - 3L^ + AL^^ + L^) - 2. (11) 

First, note that C2l+4(-^)/2 is the number of self-avoiding walks (of length 2L + 4, crossing the 
L X L square) in which the first defect, that is, the first backward, step, is a south step. We focus 
on such walks, and study four distinct cases. The first three cases count walks having two south 
steps, and the last case counts walks having a south step and a west step (Figure IH}. 
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No 



Wi 



Si 



Eo 



Ni 



Figure 8: Four types of self- avoiding walks with two defects. 



1. The walk w contains two adjacent south steps, 5i and 82- They are necessarily preceded by 
an east step Ei , and followed by another east step E2 ■ The walk has L + 4 vertical steps and 
L horizontal steps. Mark Ei, and delete 51,52, E2 in order to obtain a walk w' with L + 2 
north steps and L — 1 east steps, one of which is marked. In w' , the marked step cannot be 
at level 0, 1, L + 1 or L + 2. Using the same ingredients as above, we obtain the number of 
such walks as 

^J2L + 1\ J2L + V 



1 



L-l 



L-2 



2. The walk contains a sequence E1S1E2S2E3. Again, w has L+4 vertical steps and L horizontal 
steps. Mark Ei, and delete Si,E2,S2 and £"3 in order to obtain a walk with L + 2 north 
steps and L — 2 east steps, one of which is marked. In w' , the marked step cannot be at level 
0,l,-L-|-lorL-|-2. The number of such walks is 



(L-2) 



2L 



L 



2L 



3. The walk contains a sequence E1S1E2, and, further away, another sequence E3S2E4, disjoint 
from the first one. Again, w has L + 4 vertical steps and L horizontal steps. Mark the steps 
El and E-^, delete Si,E2,S2 and £'4 in order to obtain a walk with L + 2 north steps and 
L — 2 east steps, two of which are marked. Note that, in w' , the first marked step cannot lie 
at level 0,L + 1 or L + 2, while the second marked step cannot lie at level 0,1 or L + 2. Using 
the same ingredients as above, combined with the inclusion-exclusion principle, we find the 
number of such walks as 



'L-2 



2L 
L-2 



(L-3) 



2L 



L 



2L 
L - 4 



2L 



L 
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L-2\ f 2L \ , J 2L \ f 2L 
L-2 -^<^-"U-3 +U-4 



4. The walk w contains a sequence E1S1E2, and, further away, a sequence N1W1N2 (with obvious 
notations). It thus contains L + 2 vertical steps and L + 2 horizontal ones. Mark the steps 
El and A'^i , delete Si,E2, Wi and A''2 in order to obtain a walk w' with L north steps and L 
east steps, in which one step of each type is marked in such a way that the east marked step 
comes before the north marked step. In w' , the two marked steps cannot be consecutive (or w 
would not be self-avoiding), the east marked step cannot lie at level 0, and the north marked 
step cannot lie at abscissa L. Again, the inclusion-exclusion principle applies and gives the 
number of such walks as 



1 



2L 



{2L - 1) 



2L - 2 
L- 1 



, 2L \ f2L - r 



2L 



1 



1. 



Putting together the four partial results we have obtained gives (|TT|) . Note that the number 
C2L+4:{L) has the asymptotic expansion 

L24L f 49 2913 92971 
2- — + 



\ 4L 64L2 512L3 



The above argument suggests that it is very likely that, for every fixed K, the sequence C2l+2k{L), 
for L > 0, is polynomially recursive |161 1171 Ch. 6]. 

While it would probably be possible to find the number of possible paths of length 2L + 6, the 
number of special cases that must be treated would become onerous. We have therefore resorted 
to a numerical study for walks of length 2L + 2K, K > 2, based on our enumerations. For K = 3 
we found 

L^4L ^4 49 1931 ±1 



/L^ \3 6L 64L2 
while the corresponding result for K = 4 is 

L^A^ (2 11 



Lvr V3 4L 



We can give a heuristic argument for the general form of the leading term in the asymptotic 
expansion of the number of walks of length 2L+2i^ which gives as the leading order term —7= ^ J. . 

Here the first factor is given by the number of ways of choosing the backbone, { and the 
second is given by the number of ways of placing K defects (or backward steps) on a path of length 
2L, which is just (2L)^. The defects are indistinguishable, introducing the factor K\. 

This argument can be refined into a proof, for K = o(L^/^) by following the steps, mutatis mutandis 
in the proof of a similar result given in jHl . 
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7 Numerical analysis 



It has been proved El that hniL^oo = A exists. From this it is hkely that Rl = 

Clj^i^l+i/Cl^l ~ though this has not been proved. Accepting this, the generating function 
= "^Zl ^lx^ will have radius of convergence Xc = 1/A^, which we can estimate accurately using 
differential approximants [2]. In this way we estimate that for the crossing problem Xc = 0.32858(5), 
for the transverse problem Xc = 0.3282(6) and for the cow-patch problem Xc = 0.328574(2). It is 
reassuring to see, from our numerical studies, that A appears to be the same for the three problems, 
as proved above, and we estimate that A = 1.744550(5). 

We now speculate on the sub-dominant terms. For SAW on an infinite lattice, it is widely accepted 
that Cn ~ const. ix^n^ where c„ is the number of n step SAW equivalent up to a translation. 

It seems at least a plausible speculation that, for SAW crossing an L x L lattice, the number going 
from (0,0) to {L,L) is given by ~ AX^^^^^L"'. We have investigated this possibility numerically, 
and found it to be supported by the data, to some extent. 



We fitted the data to the assumed form, fixing the value of A at our best estimate, 1.744550. 
This then leaves two unknown parameters b and a. For cow-patch walks we find b ~ 0.8558 and 
a ~ —0.500. This suggests asymptotic behaviour ApX^ '^^■^^^^^ /^/L, and we estimate Ap ~ 0.52. 
For transverse walks and walks crossing a square b is quite small, most likely zero. A value of 6 = 
would imply the absence of a term 0(A^^), or possibly the presence of a term O(logL), or some 
power of a logarithm. We have investigated the latter possibility by including a logarithmic factor, 
and found that the data does not support the presence of such a term for either class of walk. Of 
course, we cannot rule out some small power of a logarithm, but this seems less likely than the 
absence of a term 0{X^^). 

We next investigated the possibility that the subdominant term is 0(L"). A simple ratio analysis 
P] then led to the estimates a = —0.7 for walks crossing a square, and a = 1.0 for transverse walks. 
If our assumed form is correct, we expect these estimates to be accurate to within 10-15%. We also 
studied the sequence whose terms are given by the quotient T^/Cl. This has the advantage that 
the A dependence cancels, and so our result is independent of any uncertainty in the value of A. We 
find that Tl/Cl ~ const. L^-'^ This is in agreement with the estimates of a found separately, for 
the two series. Thus we very tentatively speculate that Cl ~ 8A^ /L*^'"^ and Tl ~ 9A^ L, where 
the amplitude estimates follow by the simple expedient of fitting the assumed L dependent form 
to the data, term-by-term, and extrapolating the resulting sequence of amplitude estimates. Given 
the sensitivity of the amplitudes to both A and a, we do not feel confident quoting an uncertainty 
for the amplitudes. 



Whittington and Guttmann [T^ and later Burkhardt and Guim |2] studied the behaviour of the 
mean number of steps in a path on an L x L lattice 
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as well as the fluctuations of this quantity 



which is a kind of heat capacity. As discussed above, a phase transition takes place as one varies the 
fugacity x associated with the walk length. At a critical value x^.^ the average walk length of a path 
on an L X L lattice changes from to 0(L^). In JHj the critical fugacity was proved to satisfy 

< Xc < /U//, where ixu is the growth constant for Hamiltonian SAW on the square lattice, and 
on the basis of numerical studies conjectured to be Xc = 1/^ exactly. In ^31 the conjecture was 
proved. Here we also study the behaviour at x = Xc and find that (n(x,-L)) = 0(L^/'') where the 
numerical evidence is consistent with v = 3/4. Similar conclusions were reached earlier in [2]. For 
any given value of L the fluctuation V{x^ V) is observed to have a single maximum located at Xc{L) 
(see top left panel of Figure IHl). We study in detail the behaviour of V{x,L), which we expect to 
obey a standard finite-size scaling Ansatz 

V{x, L) ~ L2/^T/((x - (14) 

(which is equivalent to © of P]) where V{y) is a scaling function. From this it follows that the 
position and the height of the peak in V{x,L) scale as Xc{L) — Xc ~ L~^^'^ and VmaxiL) ~ L^/'^. 



Table 3: The mean-length of walks crossing an L x L square at the critical fugacity x 



the 



position, Xc{L) 



Xr 



and height, Vraax{L), of the peak in the fiuctuations V^(x, L). 



L {n{xc,L)) 



Xc{L) 



Xc 



1 2 

2 4.1230827138 

3 6.3491078353 

4 8.6519365910 

5 11.0129773423 

6 13.4187561852 

7 15.8593480600 

8 18.3273545355 

9 20.8171976528 

10 23.3246243077 

11 25.8463556412 

12 28.3798369044 

13 30.9230572826 

14 33.4744187854 

15 36.0326398605 

16 38.5966838209 

17 41.1657051788 



0.9370217352 
0.5554687338 
0.3963960508 
0.3016714640 
0.2403448999 
0.1979673072 
0.1671981710 
0.1439801106 
0.1259158112 
0.1115091953 
0.0997832765 
0.0900753740 
0.0819213689 
0.0749872153 
0.0690267737 
0.0638549420 



2.5358983849 

6.2850743202 

12.5671289312 

21.6246676036 

33.7507328831 

49.2268220069 

68.3294309970 

91.3288825240 

118.4887185709 

150.0657089122 

186.3101460060 

227.4662469752 

273.7725788463 

325.4624696518 

382.7643901657 

445.9023015941 



In table 01 we have listed the numerical values of the mean-length at Xc and the position and height 
of the maximum of the fluctuations. We analyse this data by forming the associated generating 
functions, N{z) = "^lin^x, L))z^ etc., and using differential approximants. Given the expected 
asymptotic behaviour of these quantities the generating functions should have a singularity at 
Zc = 1 with critical exponents — 1/z/ — 1 (average length at Xc), l/v — 1 (position of the peak), and 
— 2/zy — 1 (height of the peak). In table |1] we list the results from an analysis of the generating 
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Table 4: Estimates for z^. and the critical exponents obtained from second order differential approx- 
imants to the generating functions in Table 01 K is the degree of the inhomogeneous polynomial 
of the differential approximant. 



K 


(n(xc, 




XciL) 






{L) 






-Xjv - 1 


Zc 


l/v - 1 


Zc 


-2/v - 1 





0.9999823(13) 


-2.32985(17) 


1.00017(11) 


0.3147(96) 


0.999998(20) 


-3.6620(21) 


1 


0.999983(10) 


-2.3299(12) 


1.000114(23) 


0.3196(16) 


0.9999900(41) 


-3.66134(34) 


2 


0.9999818(79) 


-2.32973(99) 


1.000124(15) 


0.3185(16) 


0.999982(10) 


-3.6606(10) 


3 


0.9999789(88) 


-2.3293(11) 


1.00013(10) 


0.3183(81) 


0.999975(17) 


-3.6598(18) 


4 


0.9999773(76) 


-2.32915(93) 


1.000084(45) 


0.3215(47) 


0.999979(11) 


-3.6603(11) 


5 


0.9999786(70) 


-2.32930(80) 


1.000136(75) 


0.3171(61) 


0.9999850(69) 


-3.66081(65) 



functions using second order differential approximants. The estimates for the exponents are not 
very accurate (which is not surprising given the short length of the series) but are fully consistent 
with V = 3/4. 

Finally, in Figure El we perform a more detailed analysis to confirm the conjectured scaling form 
for y(x, L). In the top left panel we have simply plotted V^(x, L) as a function of the fugacity x to 
confirm the single peak behaviour. In the top right panel we have plotted Xc{L) and T^nax vs. L 
in a log-log plot, thus confirming that these quantities grows as a power-law with L (the straight 
lines, drawn as a guide to the eye, have slopes —l/v = —4/3 and 2/z/ = 8/3, respectively). In 
the bottom panels we check numerically the scaling Ansatz for V{x,L). In the left panel we plot 
V{x, L)/L^/^ vs. the scaling variable (x — Xc)L^^^ obtaining a reasonable scaling collapse. A better 
idea of the quality of the scaling collapse can be gauged from the plot in the bottom right panel. 
Here we plot the difference between consecutive scaling plots from the left panel. More precisely 
we plot D{x,L) = V{x, L)/L^/^ — V{x' ,L — 1)/{L — Vf/'^ vs. {x — Xc)L'^^^, where x' is chosen so 
that the scaled variables coincide, e.g., (x — Xc)L^^^ = {x' — Xc){L — 1)^/^. 



8 Walks crossing the square and hitting the centre 

In J21 Knuth also considered the problem of self-avoiding walks crossing the square and passing 
through the centre (L/2,L/2) of the grid (with L being even). Denote the number of such walks 
by c(L). Then a straightforward variant of the method of proof used in Section can be applied to 
prove that 

hm c(L)i/^' = A^. 

L— >oo 

Knuth used Monte Carlo simulations to estimate the fraction of paths hitting the centre point and 
found for L = 10 that 81 it 10 percent of all paths do hit the centre. He then went on to say that 
"perhaps nobody will ever know the true answer." Naturally we cannot let Knuth's challenge go 
unanswered. It is very simple to modify the transfer-matrix algorithm to ensure that all paths pass 
through a given vertex. We just make sure that when we do the updating at the given vertex the in- 
put state '00' (no occupied incoming edges) has only one output state '12', while the output '00' (no 
outgoing occupied edges) is disallowed at this vertex. We can thus answer Knuth's query and state 
for ah to know that for L = 10 a fraction 1243982213040307428318660/1568758030464750013214100 = 
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Figure 9: The fluctuations V{x,L) as a function of the fugacity x (top left panel). Xc{L) — Xc and 
Vmax vs. L (top right panel). V{x,L)/L^/^ (bottom left panel) and D{x,L) (bottom right panel) 
vs. the scaling variable {x — Xc)L^^^- 
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Table 5: The total number of walks crossing an L x L square and passing through the centre 
(L/2,L/2), c(L) and the ratio c{L)/C{L). 



L 


c(L) 


c{L)/C{L) 


2 


10 


0.833333 . . . 


4 


7056 


0.828947 . . . 


6 


462755440 


0.803701 . . . 


8 


2593165016903538 


0.793842... 


10 


1243982213040307428318660 


0.792972... 


12 


51166088445891978924432033203830714 


0.792927 . . . 


14 


180349587397776823066172713933745722978533730900 


0.792920 . . . 


16 


54508896286415931462305055600895616388822171335171594099162334 


0.792909 . . . 


18 


1413040380714086952244299343879218154884335669707058802937825791571640010167156 


0.792901 . . . 



0.792972. . . of all paths pass through the centre. In Table El we have listed the number of paths 
passing through the centre for L < 18. 

The fact that C{L)/c{L) appears to be going to a constant implies that not only is the asymptoti- 
cally dominant behaviour of both C{L) and c(L) the same, but so must the sub-dominant behaviour. 
We note the useful mnemonic that the ratio appears close to y^vr/S = 0.79266..., though we have 
no idea how to prove or disprove that this is the correct value. 



9 Hamiltonian walks 

Hamiltonian walks can only exist on 2L x 2L lattices. For lattices with an odd number of edges, one 
site must be missed. A Hamiltonian walk is of length 4L(L + 1) on a 2L x 2L lattice. The number 
of such walks grows as r^^ , where we find r ~ 1.472 based on exact enumeration up to 17 x 17 
lattices. In |l()j Jacobsen and Kondev gave a field-theoretical estimate of the growth constant for 
Hamiltonian SAW on the square lattice as 1.472801 it 0.00001. These were walks confined to a 
square geometry, but not restricted as to starting and end-points as are those we consider here. 
Nevertheless, it seems likely that we are estimating the same quantity, so our results can be seen 
as providing support for the view that the field theory is estimating precisely the same quantity as 
our enumerations. That is to say, this appears to be precisely the same as the corresponding result 
for Hamiltonian walks on an L x L lattice, in the large L limit. These estimates are about 20% less 
than A, the growth constant for all paths. In ^ it is proved that 2^/^ < t < 12^/^. Numerically 
this evaluates to 1.260 < r < 1.861. 

We can improve on these bounds as follows: we define cow-patch walks to be Hamiltonian if every 
vertex of the square not belonging to the border of the square belongs to one of the SAWs of 
the cow-patch. Then the upper bounds given above translate verbatim into upper bounds for r, 
while lower bounds are given by Hamiltonian traversing paths and eqn. ©. In this way we find 
1.429 < T < 1.530. As we have shown above that 1.6284 < A, this proves that r < A. 

The number of Hamiltonian paths Hi for L even, and paths that visit all but one site, for L odd, 
are given in Table El The number of Hamiltonian cow-patch paths HP^ for L even, and cow-patch 
paths that visit all but one site, for L odd, are given in Table [3 The number of Hamiltonian 
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Table 6: The number of Hamiltonian paths. 



L Hl 

1 2 

2 2 

3 32 

4 104 

5 10180 

6 111712 

7 67590888 

8 2688307514 

9 9628765945000 

10 1445778936756068 

11 29725924602729604016 

12 17337631013706758184626 

13 1998903003325610328086958408 

14 4628650743368437273677525554148 

15 2937440223891635053435045277805847436 

16 27478778338807945303765092195103685118924 

17 94555056448262478314997568263027383699860223148 



Table 7: The number of Hamiltonian cow-patch paths. 



HPl 



2 6 

3 81 

4 2420 

5 158487 

6 22668546 

7 7067228903 

8 4796951277784 

9 7083189530689311 

10 22740544515287098346 

11 158673902903632923216807 

12 2405521769596577026409223804 

13 79215226453280152797069512845071 

14 5665275864000731097175367200188234758 

15 879791999732650875090633720304683597787867 

16 296640712696590626976673730832416228749213171388 

17 217134088450048497810206709994144694071029172119163041 

18 345011492148033546292595301223727273934239259467419472922686 



transverse paths HTl for L even, and transverse paths that visit all but one site, for L odd, are 
given in Table |H1 



E-mail or WWW retrieval of series 



The series for the problems studied in this paper are available by request from I.Jensen@ms.unimelb.edu.au 
or via the world wide web http : / / www . ms . unimelb .edu.au/ ^iwan /| by following the relevant links. 
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Table 8: The number of Hamiltonian traversing paths. 



r 

L 


tllL 


1 


2 


Z 


Q 
O 


3 


34 


4 


650 


5 


12014 


6 


1 m ciACin 
iUi64yz 


7 


oo /6iyy4 


8 


32647369000 


9 


12227920752840 


10 


22181389298814376 


11 


38166266554504010420 


12 


323646210116765453608746 


13 


2574827340090912815899810042 


14 


102299512403818451392332665527950 


15 


3778748215131699995997836850757543682 


16 


704314728645701361948084580318587261484806 


17 


121135616205759617794904559766506890558675949856 


18 


106005756542854454380006180528618254764945283647525384 
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